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Abstract 

Let us denote V, the finite dimensional vector spaces of functions of the form 
i[){x) = p n (x) + f(x)p m (x) where p n (x) and p m (x) are arbitrary polynomials of 
degree at most n and m in the variable x while fix) represents a fixed func- 
tion of x. Conditions on m, n and f(x) are found such that families of linear 
differential operators exist which preserve V. A special emphasis is accorded to 
the cases where the set of differential operators represents the envelopping al- 
gebra of some abstract algebra. These operators can be transformed into linear 
matrix valued differential operators. In the second part, such types of opera- 
tors are considered and a connection is established between their solutions and 
series of polynomials- valued vectors obeying three terms recurence relations. 
When the operator is quasi exactly solvable, it possesses a finite dimensional 
invariant vector space. We study how this property leads to the truncation of 
the polynomials series. 



1 Introduction 



Quasi Exactly Solvable (QES) operators are characterized by linear differential opera- 
tors which preserve a finite-dimensional vector space V of smooth functions [1]. In the 
case of operators of one real variable the underlying vector space is often of the form 
V = V n where V n denotes the vector space of polynomials of degree at most n in the 
variable x. In [2] it is shown that the linear operators preserving V n are generated by 
three basic operators j_,j ,j+ (see Eq.(17) below) which realize the algebra s/(2,R). 
More general QES operators can then be constructed by considering the elements of 
the enveloping algebra of these generators, performing a change of variable and/or 
conjugating the j's with an invertible function, say g{x). The effective invariant space 
is then the set of functions of the form g(x)V n . 

In this paper we consider a more general situation. Let m, n be two positive 
integers and let f(x) be a sufficiently derivable function in a domain of the real line. 
Let V = V n + f(x)V m be the vector space of functions of the form p(x) + f(x)q(x) 
where p(x) e V n , q(x) G V m . 

We want to address the following questions : What are the differential operators 
which preserve V ? and for which choice of m,n,f(x) do these operators posses a 
relation with the enveloping algebra of some Lie (or "deformed" Lie) algebra ?. This 
question generalizes the cases of monomials adressed in [3] and more recently in [4] 
At the moment the question is, to our knowledge, not solved in its generality but we 
present a few non trivial solutions in the next section. 

For all the cases we obtained, the final problem of finding the algebraic modes 
of the scalar equations leads to 2x2-matrix QES operators preserving the space of 
couples of polynomials with suitable degree. On the other hand, a few example of 
systems of QES equations are known and were studied in the past. 

One problem that occur in computing the algebraic eigenvalues is to diagonalize 
the Hamiltonian in a base of the its invariant finite dimensional subspace. In the case 
of scalar QES equations, this problem can be simplified by a technique presented in 
[12]. A few years ago, Bender and Dunne [12] have pointed out that for a class of 
Schroedinger eigenvalue problems, the formal solution constructed for generic values 
of the spectral parameter E is the generating function of a set of orthogonal poly- 
nomials {P n (E)}. This property is related to the fact that the equation leads to a 
three term recursion relation on these polynomials. If the coupling constants of the 
Schroedinger operators are choosen in such a way that it becomes a quasi exactly 
solvable operator [13, 14], then the polynomials obey a remarkable property, namely 
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for n larger than a fixed integer N the polynomials P n {E) factorize in the form 
P n (E) = P^(E)P n _^(E) where Pjy is a common factor. The corresponding algebraic 
eigenvalues associated with the "quasi-exact" property are the roots of Pn{E). The 
results of [12] was generalized and considered in different frameworks in a series of 
papers [15, 16, 17, 18, 19, 20]. 

However, to our knowledge, the approach of [12] has not been applied to systems 
of quasi exactly solvable equations. Apart from their mathematical interest, such 
systems appear in a few context, e.g. in the construction of the doubly periodic 
solutions of the Lame equation [21] (see e.g. [8]) and in the stability analysis of 
classical solutions available is some low dimensional field theories [22, 23, 24],. 

It is the purpose of this letter to adapt the ideas of Bender and Dunne to systems 
of coupled Schrodinger equations and to point out that families of polynomials also 
appear in this context. The eventuality of the orthogonality of these polynomials is 
still an open question. 

The paper is organized as follows. In Sect. 2 we review two examples of well know 
coupled systems of QES equations, then we discuss the cases of operators prevers- 
ing vector spaces of the form p(x) + f(x)q(x). The constuction of solutions of the 
eigenvalue equation Hip = Eip in a form of polynomials series is adapted in Sect. 3 
for the case of coupled channels. It is shown that, along with [12] the ^-dependant 
coefficients of the series obey three terms recurence relations. Special emphasis is set 
on the relation between the finite dimensional invariant space of H and the truncation 
of the series. 

2 Systems of QES equations 

In this section, we review the most know 2x2- matrix QES equations and we study a 
class of scalar QES equations which lead, after a suitable algebra, to matrix operators 
acting on spaces of polynomials. 

2.1 Polynomial potential 

In this subsection, we study the Hamiltonian 

tf(y) = -^I 2 + M 6 (y) (1) 

where M$(y) is a 2 x 2 hermitian matrix of the form 

M 6 (y) = {4ply 6 + 8 Pl p 2 y' l + {Ap 2 1 -8mp 2 + 2{l-2e)p 2 )y 2 }lL 2 
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+ (8p 2 2/ 2 + 4pi)cr 3 - 8mp 2 K, c r i (2) 

where <7i,cr 3 are the Pauli matrices, P2,Pi,^o are free real parameters and m is an 
integer. 

It is known [9, 10] that after the standard "gauge transformation" of H(y) with a 
factor 

0(y)=Z/ e exp-{|y 4 + Pl y 2 } (3) 

and the change of variable x = y 2 , the new operator H(x) 

H{x) = ( j ) -\x)H{y) < j>{x)\ y= ^ (4) 

can be set in a form suitable for acting on polynomials. However, the invariant space 
is revealed only after the supplementary transformation [9] 

H(x) = P- 1 H{x)P , P = I 2 + ^«9 x (a! + i<7 2 ) (5) 

is performed. After a calculation, we find 

H(x) = -(Axd 2 x + 24)1 + 8pm/tofj 3 d x (6) 

+8pdiag( J + (m - 2), J + {m)) + 4/c (l + 2mpK 2 )a + d 2 x - 8mpK a_ (7) 

which manifestly preserves the space of couples of polynomials of the form 
{Pm-2{x) , q m {x)Y ■ Here we choose e = and p\ — for simplicity. If the parameter 
e is choosen as an arbitrary real number, then the initial potential M§ acquires a 
supplementary term of the form e(e — l)/y 2 - 

2.2 Lame type potential. 

As a second example, we consider the family of operators 



Ak 2 sn 2 + 5(1 + k 2 ) /2 2#A;cn dn 

20jfecndn Ck 2 sn 2 -5{l + k 2 )/2 



(8) 



where A, C, 5, 9 are constants while sn, cn, dn respectively abbreviate the Jacobi ellip- 
tic functions of argument z and modulus k [21] 

sn(z, k) , cn(z, k) , dn(z, fc) . (9) 

These functions are periodic with period AK(k),AK(k),2K(k) respectively (K(k) is 
the complete elliptic integral of the first type). The above hamiltonian is therefore 
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to be considered on the Hilbert space of periodic functions on [0,4K(k)}. For com- 
pleteness, we mention the properties of the Jacobi functions which are needed in the 
calculations 

cn 2 + sn 2 = 1 , dn 2 + A; 2 sn 2 = 1 (10) 

— sn = cn dn , — cn = — sn dn , — dn = — A; 2 sn cn (11) 
dz dz dz 

The relevant change of variable which eliminates the transcendental functions sn, 
cn, dn from (8) in favor of algebraic expressions is (for k fixed) 

x = sn 2 (z,k) (12) 

In particular the second derivative term in (8) becomes 

= 4x(l - x)(l - k 2 x)^- + 2(3k 2 x 2 - 2(1 + k 2 )x + l)-f- (13) 
dz 2 d 2 x dx 

Several possibilities of extracting prefactors then lead to equivalent forms of (8), say 

H(x), which are matrix operators build with the derivative d/dx and polynomial 

coefficients in x. The requirement that H(x) preserves a finite dimensional vector 

space leads to two possible sets of values for A, C, 9 (see [9] for details). Here we will 

discuss the set determined by A = Am 2 + 6m + 3 — 5 

C = 4m 2 + 6m + 3 + 5 
0= i[(4 m + 3)2_ 5 2]i 

The parameter S remains free, and also k which fixes the period of the potential. The 
results can be generalized easily to the other set [9]. 

Corresponding to these values, four invariant vector spaces are available. We will 
study only one of them, refering, again, to [9] for the three others. 

In order to present the change of variable we conveniently define 

R= Am + 3 - 6 , (14) 
4m + 3 + <5' K J 

We have then 



v = I J cn ° dn II ; 7 11 — 1 . * 2 = ^ 
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After performing the change of variable and the change of function, the compo- 
nents of the operator (8) take the following form 

1 5 
H u = -Axk 2 {D + m + -){D-m) + {k 2 + l){AD 2 + -)-2{l + 2D)d , 

H 12 = Axn{k 2 + l)(D-m) + k(-8D + 5 + Am + 1) , 
H 21 = + 4m + 3) , 

# 22 = -Axk 2 {D + m + ^)(D-m) 

+ (k 2 + 1)(AD 2 + 2D + 1 - -) - 2(1 + 2D)d , (16) 

where d = d/dx and D = xd. The operator H obviously preserves (V(m), (V(m)). 



2.3 Scalar QES operators preserving P n fP 7 



m 



In this subsection, we investigate [25] the linear operators preserving vectors spaces 
of the form P n © fP m and we construct several forms of the continuous function f(x), 
together with the associated values of the integers m, n for which the operators of 
interest are the envelopping algebra of some Lie (or deformed Lie) Algebra. 

2.3.1 Case f(x) = 

This is off course the well known case of [1, 2]. The relevant operators read 

U(n) = x{x± - n) , Jo (n) = ( x ± _ ^) , j_ = A. (17) 

and represent the three generators of s/(2,R). Most of known one-dimensional QES 
equations are build with these operators. For later convenience, we further define 
a family of equivalent realizations of sl(2, R) by means of the conjugated operators 
k e (a) = x a j e x~ a for e = +, 0, — and a is a real number. 

2.3.2 Case f(x) = x a 

The general cases of vector spaces constructed over monomials was first adressed in [3] 
and the particular subcase f(x) = x a was reconsidered recently [4]. The corresponding 
vector space was denoted in [4]; here we will reconsider this case and extend the 
discussion of the operators which leave it invariant. For later convenience, it is usefull 
to introduce more precise notations, setting V n = V(n, x) and 

= V {1) (N, s, a, x) = V(n, x) + x a V(m, x) 
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= span{l, x, x 2 , ... , x n ; x a , x a+1 , . . . , x a+m } 

= Vl l) ®vi 1] (18) 

in passing, note that the notations of [4] are n = s and m = N — s — 2. 

The vector space above is clearly constructed as the direct sum of two subspaces. 
As pointed out in [3, 4] three independent, second order differential operators can be 
constructed which preserve the vector space Writing these operators in the form 

J + — x(x- n)(x- (m + a)) 

dx dx 

Jo ={x- ) (19) 

. d 1 . d 
J ={x dx- + l ~ a) dx- 

makes it obvious that they preserve 

These operators close under the commutator into a polynomial deformation of the 
s/(2, R) algebra: 

[Jo,J±] =±J± 

[J+,J-] = a J 3 + (3 Jl + 7 J + 5 (20) 

where a, j3, 7, 5 are constants given in [4]. 

Clearly the operators (19) leave separately invariant two vector spaces v} 1 ^ and 
v} 1 ^ entering in (18). In the language of representations they act reducibly on 
However, operators can be constructed which preserve while mixing the two 
subspaces. The form of these supplementary operators is different according to the 
fact that the number a is an integer or not; we now adress these two cases separately. 

Case agR 

In order to construct the operators which mix and V^, we first define 

K = (D-n)(D-n + l)...D , D = x^— 

dx 

K' =(D-m-a)(D-m-a + l)...(D-a) (21) 

which belong to the kernals of the subvector spaces V n and x a V m of respectively. 
Notice that the products j e K and k e (a)K (with e = 0, ±) also preserve the vector 
space. For generic values of m, n these operators contain more than second derivatives 
and, as so, they were not considered in [4]. 

In order to construct the operators which mix the two vector subspaces entering 
in V, we first have to construct the operators which transform a generic element of 
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V m into an element of V n and vice-versa. In [5] it is shown that these operators are 
of the form 

q a = x « , a = 0,l,...,A (22) 
*« =nV(p+l-A)-j)(i A - a (23) 

where A = |m — n\ , p =max{m, n} 

The operators preserving V while exchanging the two subspaces can finally be 
constructed by means of 

Q a = q a x~ a K , Q a = x a q a K' , a = 0,1,..., A. (24) 

Here we assumed n < m, the case n > m is obtained by exchanging q a with q a in the 
formula above. 

It can be checked easily that Q a , transform a vector of the form p n + x a q m into a 
vector of the form q n G V n while Q a transforms the same vector into a vector of the 
form x a p m G x a V m - 

The generators constructed above are in one to one correspondance with the 2x2 
matrix generators preserving the direct sum of vector spaces V m © V n classified in [5] 
although their form is quite different (the same notation is nevertheless used). The 
commutation relations (defining a normal order) which the generators fullfill is also 
drastically different as we shall discuss now. First of all it can be easily checked that 
all products of operators Q (and separately of Q) belong to the kernal of the full 
space V^ l \ so we can write 

QaQf3 = Q a Qf3 = (25) 

which suggests that the operators Q's and the Q 7 s play the role of fermionic generators, 
in contrast to the J's which are bosonic (note that the same distinction holds in the 
case [5]). 

From now on, we assume n = m in this section (the evaluation of the commu- 
tators for generic values of m,n is straightforward but leads to even more involved 
expressions) and suppress the superflous index a on the the fermionic operators. The 
commutation relations between fermionic and bosonic generators leads to 

[Q, J-} =(2a-n- l)j-Q , [Q, J+] = (2a + n+ l)j+Q (26) 
[Q, J_] = -(2a + n + l)k-(a)Q , [Q, J+] = -(2a - n - l)k+(a)Q 

where the j± and k±(a) are defined in (17). These relations define a normal order 
but we notice that the right hand side are not linear expressions of the generators 
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choosen as basic elements. We also have 

[Q, D] = (D + a)Q , [Q, D] = (D - a)Q (27) 

This is to be contrasted with the problem studied in [5] where, for the case A = 0, the 
Q (and the Q) commute with the three bosonic generators, forming finally an sl(2)x 
sl(2) algebra. Here we see that the bosonic operators J and fermionic operators Q, Q 
do not close linearly under the commutator. The commutators involve in fact extra 
factors which can be expressed in terms of the operators j or k(a) acting on the 
appropriate subspace V m . This defines a normal order among the basic generators 
but makes the underlying algebraic stucture (if any) non linear. For completeness, 
we also mention that the anti-commutator {Q, Q} is a polynomial in J . 
Case a € N 

Let us consider the case a = k G N , with n < k and assume for definiteness 
m — k > n. Operators that preserve while exchanging some monomials of the 
subspace with some of (and vice versa) can be expressed as follows : 

fc— i 

W+ = x k l[(D-k-m + j) (28) 

3=0 

i n k—n—l 

W - = -U( D -j) II (D-k-n-i) (29) 

x j=0 i=l 

These operators are both of order k, W+ is of degree k while W- is of degree —k. 
When acting on the monomial of Eq.(18), W+ transforms the n + 1 monomials of 
into the first n + 1 monomials of and annihilates the k monomials of highest 
degrees in ■ To the contrary W_ annihilates the n + 1 monomials of V"/ 1 ^ and 
shifts the n + 1 monomials of lowest degrees of into . Operators of the same 
type performing higher jumps can be constructed in a straighforward way; they are 
characterized by a higher order and higher degrees but we will not present them here. 

The two particular cases k — n + 1 and k = 2 can be further commented. In 
the case k = n + 1 the space is just V m +n+i and the operators W+, W- can be 
rewritten as 

W + = ( 3+ (m + n + 1))" +1 , W. = (M n+1 (30) 

where j± are defined in (17). Setting k = 2 (and n = otherwise we fall on the case 
just mentionned), we see that the operators W± become second order and coincide 

(+2) (—2) 

with the operators noted T 2 , T 2 in Sect. 4 of the recent preprint [11]. With 
our notation they read 

W + = x 2 (D - (m + 2))(D - (m + 1)) , W_ = x~ 2 D(D — 3) (31) 
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A natural question which come out is to study whether the non-linear algebra 
(20) is extended in a nice way by the supplementary operators W± and their higher 
order counterparts. So far, we have not found any interesting extended structure. For 
example, for k = 2, n = , we computed : 

[W+, J+] = -2x 3 (D - (m + 2)){D - (m + 1))(D - m) 

[W+,J_] = -6xD(D- (m + 2))(D- -(m + 2)) (32) 

which just show that the commutators close within the envelopping algebra of the 
preserving operators but would need more investigation to be confirmed as an 
abstract algebraic stucture. 

We end up this section by mentionning that the two other vector spaces con- 
structed in [4] and denoted V^ 1 ^ and can in fact be related to by means 
of the following relations : 

V {a ~ l \x) = V {1) (N,s = 0,—^—,x a - 1 ) (33) 

CL X 

V^ a \x) = V (1 \N,s,-,x a ) (34) 

Off course the operators preserving them can be obtained from the operators above 
(19) after a suitable change of variable and the results above can easily be extended 
to these vector spaces. 



2.3.3 Case f(x) = ^Jp 2 (x), m = n — 1 



Here, P2(x) denotes a polynomial of degree 2 in x, we take it in the canonical form 
Pi{x) = (l—x)(l — \x). In the case m = n — 1, three basic operators can be contructed 
which preserve V in the case A = — 1; they are of the form 

d 

bi = nx + p 2 — 

ax 

S 2 = ^(nx-x—) (35) 

and obey the commutation relations of so(3). The operator corresponding toA 7^—1 
can be constructed from the three given above by a suitable affine transformation 
of the variable x.The family of operators preserving V is in this case the enveloping 
algebra of the Lie algebra of SO (3) in the realization above. Two particular cases are 
worth to be pointed out : 
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• A = -1 

Using the variable x = cos(0), the vector space V can be re-expressed is the 
form 

V = span{cos(n0), sin(n0), cos((n — 1)0), sin((n — 1)0), . . } (36) 

and the operators S a above can be expressed in terms of trigonometric functions. 
Exemples of QES equations of this type were studied in [6] in relation with spin 
systems. 

• A = k 2 Using the variable x = sa(z,k), (with sa(z,k) denoting the Jacobi 
elliptic function of modulus k, < k 2 < 1), and considering the Lame equation : 

-p^ + N(N+l)k 2 sn 2 (z,k)*P = E^ (37) 

It in known (see e.g. [21]) that doubly periodic solutions exist if N is a semi 
integer. If N — (2n + l)/2 these solutions are of the form 



ip(z) = \Jcn(z, k) + dn(z, k)(p n (x) + cn(z, k)dn(z, k)p n -i(x)) (38) 

where cn(z, k), dn(z, k) denote the other Jacobi elliptic functions. The second 
factor of this expression is exactly an element of the vector space under con- 
sideration. The relations between the doubly periodic solutions of the Lame 
equation and QES operators was pointed out in [8]. 



2.3.4 Case f(x) = ^(1 -x)/(l - Xx) 

For this form of f(x), the vector space V is preserved by the operators S a with 

S a = J-— S a VT^Jx~ 

V 1 — Ax 

provided m = n. 

Two cases are worth considering, in complete paralelism with Sect. 2.3.3: 

• A = —1. Using the new variable x = cos0, and using the identity tan(0/2) = 
\J(1 — x)/{\ + x) the vector space V can be reexpressed is the form 

,2n + l,, .2n+1 ,2n-l,, .2n-l, , . . 

V = span{cos( — - — 0), sm( — - — 0), cos( — - — 0), sm( — - — 0), . . .} (39) 

and exemples of QES operators having solutions in this vector space are pre- 
sented in [6]. 
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• A = k 2 . Again, in this case, the variable x = sn(z, k) is usefull and the doubly 
periodic solutions of the Lame equation (8) corresponding to N = (2n + 3)/2 
of the form [8] 

ip(z) = \Jcn(z, k) + dn(z, k)(cn(z, k)p n (x) + dn(z, k)q n (x)) (40) 
provide examples of QES solutions constructed in the space under consideration. 

2.4 Relations with matrix operators 

The scalar operators constructed in the previous subsection are equivalent to matrix 
valued differential operators acting on polynomial- valued vector. This observation is 
rather trivial but we illustrate the statement by mean of the case of Sect. 2.3.3. 

Let us now reconsider the spaces of the form V = p n + f{x)q m invariant through 
the symmetries Si, S 2 and S3 introduced above. A closed link can be obviously made 
between the later operators with matricials one which can be expressed in terms of 
the usual generators J±(n) and Jo of the s/(2,R) Lie algebra. 
We obtain iorf(x) = \J (1 — x)(l — Xx),m = n — 1 

r J_-J + {n) \ 



Si 



J_ - J+(n- 1) 



- = / -J Q {n/2)-x 2 J + {n-lY 
2 \-J (n/2) 



§= (° J-~x(D + l) 
3 [ J_ 



3 Recurence relations 

In this section, we will adapt the formulation of the QES solution in terms of recurence 
relations [12] to the case of matrix operator. Let us assume that a 2 * 2 matrix hamil- 
tonian H preserves a finite dimensional vector space. Looking for solutions of the form 

m = Mv) E Pn(E)y n , P n (E) ee ( Pn{E) ) (41) 

n=0 V ^ > I 
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We will see that the eigenvalue equation Hip = Evp leads to the a system of three 
terms recursion relations for the vectors P n (E) of the form: 

P n+1 = (El + A)P n + BP n _! 

where A and B are matrices depending on n but independant of the energy E. More- 
over, B is diagonal and posseses a zero eigenvalue for some specific value of n. We 
will see how this property, together with the form of the recurence relations lead to 
a truncation of the formal series (41). 

3.1 Polynomial potential 

The solutions of the equation 

H(x)m = Eip 

can be expressed in terms of formal series in the variable x = y 2 of the form 
i)(x) = exp(-p 2 x 2 /A) En= Pn(E)x n 

We can see after an algebra that the corresponding recurence relation has the form 






= E + B n ^ (42) 



where 




„ lib lit _L W I . . 

B n = 8p\ n (43) 



The form of C n can be obtained in a straightforward way but is irrelevant for our 
calculation. This equation can be solved recursively and it turns out that q and q\ are 
arbitrary parameters. From the structure of matrix B n ,it turns out that p m ,p m +i and 
q rn +2 can be expressed as linear combination of p m -\ and q m+ \. Therefore, considering 
the equations 

Pm-i = 0, q m +i = 

as a linear system in the free parameters qo,qi, we obtain a polynomial condition on 
the energy say P(E). Those energies that fullfill P(E) = leads to a truncated series 
and therefore correspond to the quasi-exactly solvable solutions of the initial system. 
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3.2 Lame type of potential 

After some calculus,we obtain a solution of the form 



Here,the matrices A and B are given by 



oo / i\n 

M x ) = y } tl VJE)x n (44) 



-(A; 2 + l)(4n 2 + §) -/c(-8n + 5 + 4m + 1) 
-K(5 + 4m + 3) -(A; 2 + l)(4n 2 + 2n+l-f; 

' Ak 2 {n — m — l)(n + m — \) 

4k 2 {n - m - l)(n + m + f) 



The "initial" conditions write 

6 

Pi - 1 ~ 1 



E-^(l + k 2 ) -K(5 + 4m + l 
-K(5 + 4m + 3) E - 4(1 + fc 2 )(l - 




the vector P being arbitrary. 

We immediately observe that for n — m + 1, the Matrix B vanish identically (its 
two eigenvalues are zero). Therefore all the polynomials of the form P m+ j(j=2,3,...) 
depend linearly of P m +i- As a consequence, solving the system of two equations 
P m +i = as function of the two arbitrary constants entering in P Q leads to a trunca- 
tion of the series and a corresponding set of algebraic eigenvalues. This finally leads 
to a factorization of all the polynomials P m+ j(j=3,4,...) in terms of Pj+i- 



3.3 Bose-Hubbard Model 

In the above sections, we illustrated some new aspects of the formulation of the QES 
property in terms of recursive polynomials. Here we would like to present another 
exemple where still a new feature of this formulation is revealed. It is connected to 
the Bose-Hubbard (BH) model. It can be described by the Schrodinger equation as 

[-^ + V{x)]4>{x) = E^ (45) 

where V(x) = ^ cosh 2 ( v /7:r) — (n+1) cosh( v /7x) — ^— 7f (f — 1) and | is the coefficient 
of the term aWaa term in the BH Hamiltonian of Ref. [30, 31]. 
The Eq. (45) can be written as 

d 2 1 

[— -— + (— cosh arc — M) 2 ]ip = Eip (46) 
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where, a = ^f, E = E 1 + E , E = ^-a 2 + £ + ^n(n — 2), M — 

Now to see the orthogonal polynomial associated with this model we substitute, 



ip{x) = exp ( cosh ax)4>(x) (47) 



a 2 



in the eq. (46) to obtain, 



2 1 2M 1 
0" sinh ax (j)' +[E — M 2 + ( ^ ) cosh ax\(j) = (48) 



a a 2 a a 2 ' 



The Eq.(48) further can be reduced to, 

a 2 z(z + 2) ( j ) "+[a 2 (z + l)-2z 2 -4z} ( j ) '+[E-^-M 2 + (— - i-)(^+ 1)]0 = (49) 

cr a cr 



where z = cosh ax — 1 and primes indicate derivatives with respect to z. Now this 
equation has regular singular point at z — 0, therefore we seek a solution of the form, 

4>(z) = z s f(z) (50) 

By substituting this in Eq. (49) and putting the coefficient s of the term z s ~ 1 f(z) 
equal to zero, we obtain the indicial equation for s as 2s 2 = s . This implies s can be 
either or |. The differential equation in terms of f can be written as, 

a 2 [(z + 2) 2 - 2(z + 2)] /" + [-(z + 2) 2 + (z + 2){a 2 (2s + 1) + 2} - a 2 } f 



2M ,2M 1 



+ [E-M 2 + a 2 s 2 + ( -)(z + 2)]/ = (51) 

Q: Ol Ol 



We further substitute 



in the Eq. 51 to obtain the three term recursion as 

ry 2 n 2 n 2 2M 

^-R n+2 (E) = R n (E)[E + ^- + sna 2 + n + s 2 -M 2 -—} 
4 4 a 

+ R n _ 2 (E)n(n-l)[—-\-n-2] (53) 

provided 2s 2 = s. Thus we have two sets of independent solutions: the even states 
(i.e., states with even number of nodes) for s = and the odd states for s = \. 
Note that unlike Bender-Dunne ( or most other QES) cases, s is not contained in the 
potential and this is perhaps related to the fact that for any integer value of M with 

~ 2aM - 1 
M= ^— 

cr 
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the QES solutions corresponding to both even and odd states are obtained. Also from 
Eq, 53 we observe that the even and odd polynomials, R n (E) do not mix with each 
other and hence we have two separate three-term recursion relations depending on 
whether n is odd or even. In particular, it is easily shown that three term recursion 
relations corresponding to the even and odd n cases , respectively are given by, n > 1 



op- 2M 

—P n {E) = Pn-AE) [E + a 2 (n 2 - 2n + 1 + 2ns - 2s) + 2n - 2 + s 2 - M 2 ] 

4 a 

+ 2P n _ 2 (£) („_i)(2n-3)[— -Ar-2n] (54) 

a cr 

a 2 1 2M 

—Q n (E) = Q n -i(E) [E + a 2 (n 2 -n + - + 2ns-2s) + 2n-l + s 2 -M 2 ] 

4 4 a 

+ 2Q n _ 2 (E) ( n -l)(2n-l)[— -\-2n-l] (55) 

a cr 

with Pq(E) = 1 Qo(E) = 1. These recursion relations generate a set of monic 
polynomials and forms separately the complete set orthogonal polynomials. These 
odd and even polynomials satisfied the factorization properties of the Bender-Dunne 
polynomials. It is easily seen that when M is a positive integer, exact solutions for 
first M levels are obtained . In particular, if M is odd (even) integer, then solutions 
with even number of nodes (s = 0) are obtained when the coefficient of P( n _2) (Q( n -2) 
vanishes. Similarly if M is odd (even) integer, the solution with odd number of nodes 
(s = |) are obtained when the coefficient of Q( n -2) {P{n-2) vanishes. Further for M 
even (say 2k + 2, k — 0, 1, 2 • • •), half of levels, i.e. k + 1 levels are obtained each from 
the zeros of the orthogonal polynomials Pk+i(E) and Qk+i(E) . On the other hand , 
when M is odd (say 2k + 1, k — 0, 1, 2 • • • ) then k + 1 and k levels are obtained from 
the zeros of the orthogonal polynomials Pk+i(E) and Q k (E), respectively. 



4 Concluding remarks 

The examples of operators presented above give evidences of the difficulty to classify 
the coupled-channel (or matrix) QES Schrodinger equations. The way of constructing 
the QES potential M 6 in Sect. 2 further provides a clear link between the approaches 
[26] and [28] to this mathematical problem; we hope that this note will motivate 
further investigations of it. We also presented a few new aspects of the recurence 
relations for polynomials associated to QES operators. Especially we formulated it 
for matrix operators. An open question is to show that these vector- valued polynomial 
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are orthogonal with respect to an appropriate measure. We failed to find the matrix 
counterpart of the orthogonality theorem [32] for vector-valued polynomials. 
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